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Abstract 

Context. Molecular lines in the (sub)millimeter wavelength range can provide important information about the physical and chemical 
conditions in the circumstellar envelopes around Asymptotic Giant Branch stars. 

Aims. The aim of this paper is to study the molecular composition in the circumstellar envelope around the oxygen-rich star IK Tau. 
Methods. We observed IK Tau in several (sub)millimeter bands using the APEX telescope during three observing periods. To deter- 
mine the spatial distribution of the '-C0(3 - 2) emission, mapping observations were performed. To constrain the physical conditions 
in the circumstellar envelope, multiple rotational CO emission lines were modeled using a non local thermodynamic equilibrium 
radiative transfer code. The rotational temperatures and the abundances of the other molecules were obtained assuming local thermo- 
dynamic equilibrium. 

Results. An oxygen-rich Asymptotic Giant Branch star has been surveyed in the submillimeter wavelength range. Thirty four tran- 
sitions of twelve molecular species, including maser lines, were detected. The kinetic temperature of the envelope was determined 
and the molecular abundance fractions of the molecules were estimated. The deduced molecular abundances were compared with 
observations and modeling from the literature and agree within a factor of 10, except for SO2, which is found to be almost a factor 
100 stronger than predicted by chemical models. 

Conclusions. From this study, we found that IK Tau is a good laboratory to study the conditions in circumstellar envelopes around 
oxygen-rich stars with (sub)millimeter-wavelength molecular lines. We could also expect from this study that the molecules in the 
circumstellar envelope can be explained more faithful by non-LTE analysis with lower and higher transition lines than by simple LTE 
analysis with only lower transition lines. In particular, the observed CO line profiles could be well reproduced by a simple expanding 
envelope model with a power law structure. 

Key words, asymptotic giant branch star - molecules - abundances 
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1 . Introduction The chemistry of the atmospheres and, further out, of the cir- 

cumstellar envelopes (CSEs) around AGB stars is dependent on 
^. Stars with initial masses lower than ~8 evolve to a pulsa- chemical class. They are classified either as M stars (C/O 
^ tionally unstable red giant star on the Asymptotic Giant Branch abundance ratio < 1), S stars (C/O ^ 1) or C stars (C/O > 1). The 
O . (AGB). At this stage, mass loss from the evolved central star -^^^ -^^^^^^ ^ ^^^^^ ^1^^^ absorption from 
O . produces an expanding envelope. Furt her on, carbon C is ^^^jj^^ atmosphere. M-type stellar spectra are dominated 
• • "fp I 'I T\loA^ ' |y^'"^"i"'-^^tal. || l996t by oxygen-bearing molecules, e.g., the metal oxides 
> [Fukasaku et al. 1994). r^^Q j^^q atmospheres carbon-bearing 
AGB stars are characterized by low surface temperatures, j^olecules like, a.o., CH, C2, C2H2 and HCN are detected at op- 
r. < 3000 K, high luminosities up to several 10^ U, and a very ^-^^^ -^^^^^^^ wavelenths, and in the microwave regime (e.g. 
large geometrical size up to several AU ( |Habmsl | 1996|) . In gen- Oautschv-Loidl et al. 2004). While the atmospheric abundance 
eral, these highly evolved stars are surrounded by envelopes with fractions are nowadays quite well understood in terms of ini- 
expansion velocities between 5 km^s-' and40kms-i. They have ^^^j chemical composition, which may be altered by nucleosyn- 
high mass-loss rates between 10 and 10 Mq yr . Their at- ^^etic products which are brought to the surface due to dredge- 
mospheres provide favorable thermodynamic conditions for the ^p^^ ^^e main processes determining the circumstellar chemical 
formation of simple molecules, due to the low temperatures and, abundance stratification of many molecules are still largely not 
simultaneously, high densities. Due to pulsation, molecules may understood. In the stellar photosphere, the high gas density en- 
reach a distance at which the temperature is lower than the ^^^^^ ^j^^^^^^j equilibrium (TE). Pulsation-driven shocks in the 
condensation temperature and at which the density is still high ^^^^^ ^^^^ ^egi^^ suppress TE. This region of strong shock ac- 
enough for dust grains to form. Radiation pressure drives the jj^j^y ^j^^ the locus of grain formation, resulting in the deple- 
dust away from the star Molecules surviving dust formation are ^-^^ molecules as SiO and SiS. Other molecules, as CO 
accderated due to dust-grain colhsions ( |Goldreich & Scovillej and CS, are thoug ht to be inreactive in the dust forming region 
| 19/q) . jDuari et al.lll999t) . At larger radii, the so-called outer envelope 
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is penetrated by ultraviolet interstellar photons and cosmic rays 
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resulting in a chemistry governed by photochemical and ion- 
molecule reactions. This picture on the chemical processes alter- 
ing the abundance stratification is generally accepted, but many 
details on chemical reactions rates, molecular left-overs after the 
dust formation, shock strengths inducing a fast chemistry zone 
etc. are not yet known. 

Spectroscopical studies of molecular lines in the 
(sub)millimeter range are a very useful tool for estimating 
the physical and chemical conditions in CSEs. Due to its 
proximity, the carbon-rich AGB star IRC-i-10216 has attracted 
lot of attention, resulting in the detection of more than 60 
differ e nt chemical conipound s in its CSE (e.g. iRidgwav et al.l 
Il976t ICernicharo et all l2000h . Until now, detaile d studies of 
oxyge n-rich envelopes have been rare. Recently, IZiurvs et alJ 
(|2004|) have focused on the chemical analysis of the oxygen-rich 
peculiar red supergiant VY CMa. VY CMa is, however, not a 
proto-type of an evolved oxygen-rich star A compl ex geometry 
is ded uced from Hubble Space Telescope images (ISmith et alj 
1200 Ih with a luminosit y larger than 10^ Lq and a mass-loss rat e 
of ~2 X 10^4 Mjyj. dBowers etalJ 119831: LSooka et al.l flQSSh . 
VY CMa is a spectacular object, which because of its extreme 
evolutionary state can explode as a supernova at any time. 
Interpreting the molecular emission profiles of VY CMa is 
therefore a very complex task, subject to many uncertainties. To 
enlarge our insight in the chemical structure in the envelopes of 
oxygen-rich low and intermediate mass stars, we therefore have 
started a submillimeter survey on the oxygen-rich AGB star 
IK Tau, which is t hought to be ( roughly) spherically symmetric 
dLane et al.lll987l: iMarvel 120051) . We thereby will advance the 
understanding on the final stages of stellar evolution of the 
majority of stars in galaxies as our Milky Way and their resultant 
impact on the interstellar medium and the cosmic cycle. 



1.1. IK Tau 

The Mira variable IK Tau, also known as NML Tau, is lo- 
cated at Qr20oo=3''53™28^8, 520oo=ll°24'23". It was found to 
be a n extremely cool star having large infrared {J - K) excess 
dXicolea et al. 1999) consistent with a 2000 K blackbody. IK Tau 
shows regular optical variations with an amplitude of ~ 4.5 mag. 

IK Ta u is an O-rich star of spec tral type ranging from M8. 1 
to Ml 1.2 dWing & Lockwoodll973h . Its distance was derived by 
lOlofsson et alJ (11998 ) to be 250 pc assuming a stella r tempera- 
ture o f 2000 K. The pulsation period is ~470 days dHale et al.l 
Il997h . The systemic velocity of the star is 33.7 kms"'. Mass- 
loss rate estimates range from 2 . 4 x 10 '^ Moyr"' (from the 
CO(J=1-0) line; lOlofsson et all 1 1998h to 3 x10'^ M^yr ' 



Table 1. Overview of the molecular line transitions observed 
with APEX. 



(from an analysis of multiple SiO lines; iGonzalez Delgado et al.l 
120031) . 

In the circumstellar envelo pe of I K Tau maser em ission 
from OH dBo wers et al.' '19891), H2O (iLane et al.l [19871) . and 
SiO (iBoboltz & Diamond 2005) and thermal emission of SiO, 
CO, SiS, SO, S O? a n d HCN have pr e viousl y been found 



dLindqvist et al.l 119881; iBuiarrabal etal] Il994t lOmont et al.l 
119931) . Obviously, IK Tau is a prime candidate for circumstel- 
lar chemistry studies. 

2. Observations 

The observatio ns were performed with the APExQ 12 m tele- 
scope in Chile (iGiisten et al.ll2006l) located at the 5100 m high 
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Notes. The third column lists the transition frequency, the last column 
the beam size (HPBW). For CN, only the strongest hyperfine compo- 
nent is given. 



site on Llano de Chajnantor. The data were obtained during 
observing periods in 2005 November, 2006 Apr il and August. 
The r eceivers used were the facility A PEX-2A ( Risacher et al.l 
I2OO6I) and the MPIfR FLASH receivers dHevminck et al.ll2006l) . 
Typical system noise temperatures were about 200 K - 1000 K 
at 290 GHz and 350 GHz, 1000 K at 460 GHz and 5000 K at 8 10 
GHz, respectively. The spectrometers for the observations were 
Fast Fourier Transform Spectrometers (EFTS) with 1 GHz band- 
width and the channel width for the 290-350 GHz observations 
was approximately 122.07 kHz (8192 channels), and for the 460 
GHz and 810 GHz observations 488.28 kHz (2048 channels). 
For the observations, a position-switching mode was used with 
the reference position typically 180" off-source. The antenna 
was focused on the available planets. IK Tau itself was strong 
enough to serve as a line pointing source, thus small cross scans 
in the '^CO(3-2) line were done to monitor the pointing during 
the observations. The telescope beam sizes (HPBW) at frequen- 
cies of the observed molecular lines are shown in Table [1] The 



' This publication is based on data acquired with the Atacama 
Pathfinder Experiment (APEX). APEX is a collaboration between 



the Max-Planck-Institut fur Radioastronomie, the European Southern 
Observatory, and the Onsala Space Observatory. 
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Table 2. Beam efficiencies for the different receivers. 



IK Tau 



Receiver 


Forward efficiency 


Beam efficiency 


APEX-2A 290 GHz 


0.97 


0.80 


APEX-2A 350 GHz 


0.97 


0.73 


FLASH 460 GHz 


0.95 


0.60 


FLASH 810 GHz 


0.95 


0.43 




50 

Velocity (km/s) 



Figure 1. Observed '"CO and '^CO line profiles (solid lines) to- 
gether with the expanding shell fit (dotted lines). 



antenna beam efficiencies are given in Table l2l of iGiisten et al.l 
(I2006h . 

To map the circumstellar envelope in the '^CO(3-2) line, 30 
positions distributed on a 5x6 grid in right ascension and decli- 
nation were observed. The grid spacing was 9" (half the FWHM 
beam size at 345 GHz). A raster mapping procedure was used 
along the parallel grid lines with an integration time of 15 s. 

The spectra were reduced using the CLASS program of the 
IRAM GILDASB To calculate the main-beam brightness tem- 
peratures of the lines, Tmb, the following relation was used: 



MB 



7,* If 
^ A ■ 



(1) 



Here is the measured antenna temperature, rjf is the for- 
ward efficiency and ri^ff is the antenna main-beam efficiency of 
APEX (see Table|2|. 



^ GILDAS is a collection of softwares oriented toward (sub- 
)millimeter radio astronomical applications developed by IRAM (see 
more details on http ://w ww.iram.fr/IRAMFR/GILDAS ) . 




Figure 2. Observed HCN, ^''SiO and ^"SiO line profiles (soUd 
lines) together with the expanding shell fit (dotted lines). 



3. Observational results 

Thirty four transitions from 12 molecular species including 
maser lines were detected with the APEX telescope toward IK 
Tau. The detected molecular lines are listed in Table [T] and their 
spectra are displayed in Figs.[T]to|5] 

Fig. |6] and Fig. |7] show the H2O maser lines and SiO maser 
lines observed toward IK Tau, respectively; the maser line pa- 
rameters are given in Table [3] Maser emission from H2O at 
321 GHz and 325 GHz was detected, as well as in the 7 = 7-6 
and 7 = 8-7 rotational transitions within the v = 1 and v = 3 
vibrationally excited states of -**SiO, ^"^SiO and ^^"SiO. 



3.1. Line parameters 

To get the mean brightness temperature estimates, the spectra 
were correct ed by the beam- filling factor s assuming a CO source 
size of 17" (Buiarrabal & Alcoleal[T99l . a HCN source size of 
3.85" (|Marv el 2005) an d source sizes for the other molecules of 
2.2" dLucas et ani992l) . Note that the CO size may be uncertain, 
likely underestimated, since the signal-to-noise (S/N) ratios of 
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Velocity (km/s) Velocity (km/s) Velocity (km/s) 

Figures. Observed CN, SO and CS line profiles (solid lines) together with the expanding shell fit (dotted lines). 
Table 3. Line parameters for the detected maser lines. 
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1.22 


33.4 
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0.41 
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^"SiO 
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1.40 


5.4 


33.8 



Table 4. Line parameters for each (non-masering) transition. 
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32.7 


9.01 


0.07/4.76 
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18.7 



Notes. The third column gives the upper energy level and the fourth column the line strength. The fifth and sixth column yield the peak values and 
integrated intensity values derived from the observed line profiles. In the last column, the expansion velocity is derived from the expanding shell 
fit. For the peak flux and the integrated intensity, the values derived from the main-beam temperature (Tmb) and the mean brightness temperature 
estimate (Mean Tb) are given. 



4 



H. Kim et al.: The circumstellar chemistry in the O-rich AGB IK Tau 



0,2 



0.2 



0,5 



0,5 



SiS(16-15) 





50 
Velocity (km/s 



100 



Figure 4. Observed SiS and SiO line profiles (solid lines) to- 
gether with the expanding shell fit (dotted lines). 



the profiles obtained by iBuiarrabal & Alcoleal (Il99lh are much 
smaller than those of the CO profiles presented in this paper 
The beam-filling factor is given by 



(2) 



where ds is the source size and 9b is the half-power beam width 
(HPBW) shown in Table[T] Both source and beam are assumed to 
be circular Gaussians. The mean brightness temperature estimate 
is computed by 



_ 1 

Tb - y^^MB . 




50 
Velocity (km/ s) 



100 



(3) 



Figure 5. Observed SO2 line profiles (solid lines) together with 
the expanding shell fit (dotted fines). 



Line parameters were derived with CLASS (see more details on 
http://www.iram.fr/IRAMFR/GILDAS) from fitting the spectral 
lines with expanding shell fits, from which the expansion veloc- 
ity of the envelope is obtained. The observed maser line and ther- 
mal emission line parameters are given in Table[3]and|4] includ- 
ing the envelope expansion velocity Vexp, the main beam bright- 
ness temperature Tub, the integrated area, and the parameters of 
the expanding shell fits. The expansion velocities are distributed 
from 14 kms"' to 21 kms"'. 

When the S/N ratio is high enough to warrant a considera- 
tion of the shape of the line profiles, they appear to be character- 
istic for circumstellar envelopes (for more detail see Zuckermaril 
Il987h : the '^CO lines have the parabolic shape of optically thick 
lines and the '^CO(3 - 2) line has the double-horn shape of spa- 
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Figure 6. 321 GHz and 325 GHz H2O maser emissions observed 
towards IK Tau. 
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Figure 7. ^^SiO, ^^SiO and ^''SiO maser emission detected in the 
rotational transitions 7 = 7-6 and 7 = 8 - 7 of the fundamental 
v=l and v=3 vibrational excited states observed towards IK Tau. 



tially resolved optically thin lines (see Fig. [TJ. Lines from the 
three SiO isotopologues and SiS lines have a Gaussian shape 
(see Fig.|4]i, indicating that they are partially formed in the wind 
acceleration regime wher e the stellar winds has not y et reached 
its full terminal velocity jBuiarrabal & Alcolealll99lT) . Some of 
the SO2 lines seem to show the square shape characteristic of un- 
resolved optically thin lines and some of them have the parabolic 
shape of optically thick lines (see Fig. |5]l. CS and SO lines 
seem to have the square shape of unresolved optically thin lines 
for low excitation transitions and the parabolic shape of opti- 
cally thick lines for high excitation transitions (see Fig.[3]l. HCN 
shows a global parabolic shape with a weak double-peak profile 
on the top (see Fig. |2]i. For the CN molecule fits to the spec- 
tra were done that take the hyperfine structure of the molecule 
into account. Although the S/N of the individual components is 
small, the observations are not in agreement with the optical thin 
ratio of different HFS co mponents and hint to hy perfine anoma- 
lies as already reported bv lBachiller et al.l(ll997h . 



3.2. CO Maps 

The spectra resulting from mapping the '^CO(3 - 2) transition 
in a region of 45"x54" around IK Tau are shown in Fig. [8] 
These spectra provide us with a tool to derive the source size 
as a function of radial velocity (see Fig. |9j.The envelope of IK 
Tau appears roughly spherically symmetric in '^CO(3 - 2) with 
a deconvolved extent at half-peak integrated intensity of 20". 
The physical diameter of the emission region is thus 2.1 x 10'* 
cm assuming a source distance of 250 pc. 



4. Modeling results 

4.1. Physical structure of the envelope 

CO lines are amongst the best tools to estimate the global prop- 
erties of circumstellar envelopes, since the abundance of CO is 
quite constant across the e nvelope, except for p hoto-dissociation 
effects at the outer edge (iMamon et al.lll98"8l) . The spatial dis- 
tribution of CO was found from our mapping observation to be 
spherically symmetric (see Sect. 13.21 ). A detailed multi-line non- 
LTE (non local thermodynamic equilibrium) study of CO can 
therefore be used to determine the physical properties of the en- 
velope. 

The o ne-dimensional version of the M onte Carlo code 
RATRAN dHogerheiide & van der T^l2000l) was used to sim- 
ulate the CO lines' emission. The basic idea of the Monte 
Carlo method is to split the emergent radiative energies into 
photon packages, which perform a random walk through the 
model volume. This allows the separation of local and exter- 
nal contributions of the radiation field and makes it possible 
to calculate the radiative transfer and excitation of molecular 
lines. The Monte Carlo method f or molecular line transfer has 
been described by iBernes 1 (Il979h for a spherically symmetric 
cloud with a uniform density. The code is formulated from the 
viewpoint of cells rather than photons. It shows accurate and 
fast performance even for high opacities (for more details see 
[hogerheiide & van der Tak 2000). The circumstellar envelope is 
assumed to be spherically symmetric, to be produced by a con- 
stant mass-loss rate, and to expand at a constant velocity. In the 
Monte Carlo simulation, typically 10^ model photons are fol- 
lowed throughout the envelope until they escape. The region 
is divided into discrete grid shells, each with constant proper- 
ties (density, temperature, molecular abundance, turbulent line 
width, etc.). 
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Figures. Map of '^CO J=3-2 line emission toward IK Tau. The grid spacing was 9". The main-beam brightness temperatures in 
the spectra range from 0.5 K to 3.3 K. 



For the case of a steady state, spherically symmetric outflow, 
the gas density as a function of radial distance from the center of 
the AGB star is given by 



n{r) 



M 



(4) 



where m is the mass of the typical gas particle, taken to be m ~ 
3 X 10"^"* g ram, since the gas is m ainly in molecular form in AGB 
envelopes dTevssier et alj|20()6l) . 

The kinetic temperature is assumed to vary as 



10 



16 



r(cm) 



(5) 



where Tq is the temperature at 1 x 10"' cm and Tbg represents 
the background temperature. With the radial profiles for den- 
sity and temperature given by Eq. |4] and |5j the program solves 
for the molecular excitation as a function of radius. Beside col- 
lisional excitation, radiation from the cosmic microwave back- 
ground and thermal radiation from local dust were taken into 
account. Then, the molecular emission is integrated in radial di- 
rection over the line of sight and convolved with the appropriate 
antenna beam. 

The best-fit model is found by minimizing the total x~ using 
the;^^- statistic defined as 



X 



[^mod ^obs ] 



(6) 



where /mod is the line intensity of the model and /obs is the ob- 
servation, cr is rms noise of the observed spectra, the summation 
is done over all channels of the three '^CO line transitions as 
observed for this project with APEX, i.e. J - 3-2, J - 4-3, and 
7 = 7-6. We have put more weight on the reproduction of the 
line shapes and the fitting of the lines observed with the APEX 



telescope which were calibrated in a consistent way than on the 
reproduction of the lines taken from the literature. The reduced 
for the models is given by 



X 



(7) 



where d.f is the degree of freedom being N - p, with p the 
number of adjustable parameters. Fig. [TOl shows the;^'^ contour 
plot produced by varying the mass-loss rate and the temperature 
To. In this figure, the 68 % confidence limit, i.e. the Icr level, is 
indicated. In this region, the temperature Tq ranges between 34 
to 47 K and the mass-loss rate is in the range from 4.0x10"*' to 
5.7x10-^ Mo/yr 

The best-fit model parameters are listed in Table |5j the re- 
sults of the model fits are shown in Fig.[TT] In Fig. fT2l theoretical 
model predictions for the '-CO lines with different inner radii, 
different To, and different outer radii are shown. Predictions for 
'^CO with different Tq are presented in Fig.[T3] Predictions for 
intensities at the observed offset positions were done from the 
best fit model and are consistent with the size determined from 
the observed CO maps. 

As shown in Fig. [TT] the overall line profiles are fit very 
well for the higher / transitions ('-CO(3 - 2), ^^CO{4 - 3), 
'^CO(7 - 6)). However, the model intensities of the IRAM 
'^CO(l - 0) and '^CO(2 - 1) transitions are somewhat higher 
than the observational data taken from literature, but the shapes 
fit satisfactorily. The predictions for the '^CO(2 - 1) line are still 
within the absolute uncertainty of the line, but this is not the case 
for the '-CO(l - 0) line. An obvious reason for this mismatch 
could be a problem with the outer radius value. However, our 
sensitivity analysis (see Fig.[T2land see discussion in next para- 
graphs) shows that, while lowering the outer radius value indeed 
the total integrated intensity decreases, the line shape is not well 
reproduced anymore. Since the relative uncertainty (i.e., the line 
shape) is much lower than the absolute intensity (i.e., the inte- 
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Figure 9. Contour maps of the '^CO i-3-2 line emission for IK Tau integrated over 10 km s ' centered at velocities of 15 km s 25 
kms"', 35 kms"' and 45 km s"'. One velocity interval is centered near the stellar radial velocity of 33.7 km s '.The contour values 
are 20 % to 90 % and 99 % of the maximum integrated intensity in each velocity interval, which is 7, 29, 32, and 26 Kkms ' for 
the 15, 25, 35, and 45 km s ' channel, respectively. The 50 % contour level is drawn in boldface. 




? la 4 10 _e 10 

Moss-loss rote (M^^i^^ yr ) 

Figure 10. Reduced contour maps of the model to input 
variable parameters, i.e., the mass-loss rate (M) and temper- 
ature (To), which is the kinetic temperature at a given radius 
(1 X 10'^ cm). The contour level is drawn at 0.94 indicating the 
68 % confidence level. 



grated intensity), we have put more weight on the reproduction 
of the line shapes. Moreover, we note that this is not the first 
time that a non-compatabili ty of the IRAM flu xes with other ob- 
served data is reported (e.g. lDecin etaDl2008h . The '^^CO(3 - 2) 
line clearly shows a double-horn profile and the best-fit results in 
a somewhat different Tq and a different outer radius than for the 
'^CO data. Nevertheless, the best-fit value for To derived from 



'^CO still gives a reasonable fit to the '^CO line (Fig. [T3T l. As 
shown in Fig. [13] the intensities of the profiles do not change so 
much with Tq but the lines show a flat shape on top for the lower 
temperatures (30 K and 40 K), and a double-horn shape at higher 
temperatures. 

As shown in Fig. \T2\ the line shapes and intensities for all 
transitions are not much influenced by the inner radius variations 
since the '"CO emission contributing dominantly to the spectra 
arises from regions further out in the envelope. The outer radius 
variations mainly aff'ect the '^CO(l -0) line, which is formed 
further out in the envelope than the other '^CO transitions. 

4.2. Chemical abundance structure 

As explained in the introduction, the density distribution of each 
molecule is different, depending on the chemical processes par- 
taking in the envelope. The fractional abundance of a species A 
is usually specified as 



Mr) = 



(8) 



where «h,('") is the number density of H2 and«A('") is the number 
density of species A. 

A first order assessment of the molecular abundance frac- 
tions can be obtained assuming that the envelope structure is 
in local thermodynamic equilibrium. Assuming a spherically 
symmetric envelope, the fractional abundance for an optically 
thin rotational line {J — » 7 - 1) of a linear rotor is given by 
lOlofsson etai] ( ll991h . 
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Figure 11. Best-fit model spectra for the different CO transitions to the available set of data at offset (0",0") for IK Tau. 
Parameters used to calculate the best-fit theoretical predictions are given in Table |5] The '-CO(l - 0) and '-CO(2 - 1) data are 
from .Tevs sier et al. (2006), where 'pv' means the IRAM 30 m telescope in Pico Veleta. The model fits are shown by dotted lines. 





^ 1 1 nt i 1 1 1 I ^Vt- I 

6 P ''cO(2-l)-pv 

Tj 
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Figure 12. Model fits with different input parameters to the set of data at offset (0", 0") for IK Tau. Left column shows variations 
in the inner radius: 0.1 x 10'"^ cm (dash-dotted line), 1 x 10'"* cm (dashed line), 10 x 10'"* cm (dotted line). Middle column shows 
variation in Tq: 30 K (dash-dotted line), 40 K (dashed line), 50 K (dotted line). Right column shows variations in the outer radius; 
5.3 X 10"" cm (dash-dotted Une), 6.3 x 10"" cm (dashed line), 7.3 x 10'^ cm (dotted Une) . 
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Table 5. Parameters for the best-fit model to the observed '^CO and '^CO line profiles. 





Roul 


Mass-loss late 


fco 


V 

* exp 


To 


a 




(lO''* cm) 


(lO''* cm) 


(Moyr-') 


(10-*) 


(kms-') 


(K) 




(K) 


'-CO 1 


630 


4.7x10-" 


3 


18 


40 


0.8 


2.7 


'^CO 1 


700 


4.7x10-** 


0.35 


18 


50 


0.8 


2.7 



Notes. Rj represents the inner radius of the envelope, R^ui the outer radius and fco=[CO/H2] the CO abundance relative to H2. The last three 
columns give the values for the parameters in Eq.|5] The expansion velocity is derived from the '^CO(3-2) observations. 




Figure 13. Model fits of the '^CO(3-2) transition with different values for Tq, as indicated in the upper left corner of each panel. 
The model fits are shown by dotted lines. 



V- BD T 
3.8x10 T^,—^ — 

/i(Z?oXlOV(^+l), 
• exp( — ) 



J^' exp(-4 In 2x^)dx 



(9) 



where Tmb is the main-beam brightness temperature, T^x is 
the excitation temperature {- Tiot and equal to the kinetic tem- 
perature under the LTE assumption), fio is the dipole moment 
in Debye, Bo is the rotational constant in GHz, Vgxp is the gas 
expansion velocity of the CSE in kms"', B is the beam size in 
arcseconds, D is the distance to the source in pc, M is the mass- 
loss rate in Moyr ', and x, (=0) and are the inner and outer 
radius of the CSE, respectively, measured in units of B. It has for 
simplicity been assumed that f ^ is constant from Xi to x^ and zero 
elsewhere. For CN, the relative strenghts of the different hyper- 
fine components where taken into account. If the line is optically 
thick, the value of /v estimated by the above formula is only a 
lower limit. 

The SO2 abun dance with respect t o H2 is estimated using the 
equation given bv lMorris et al. I (ll987l) : 



/SO2 - 



2x1013 M^lSv 

1 1 
exp[-1.44£u/rex] exp(-41n2x2)£/jc 



(10) 



where QiT^x) i s the molecular partit ion function (^ l.lST^l^, 
for more detail see lOmont et al.l (Il993h ). fiu is the energy of the 



upper state of the transition, S is the line strength, and v is the 
frequency of the transition. 

A mass loss ra t e of M - 4.7 x 10"'' Moyr"' (see Sect. 14. II 
and iTevssier et al.l l2006l) was adopted to calculate the abun- 
dances. Since the outer radius of the molecular emitting re- 
gion can be quite uncertain for molecules for which no observa- 
tional maps exist, two different outer radii will be used for these 
molecules ('case A' and 'case B'). For SiO, the value for the 
outer radius was taken to be 2 x 10'^ cm (case A) and 5 x 10'^ cm 
(case B), for the other mole cules Ix 10'^ cm (case A ) andSxlO'^ 
cm (case B) was assumed jBuiarrabal et al.l 19941) . For all lines 
from this work, we adopted expansion velocities from Table |4] 
For lines taken from the literature (see Table [T]), an expanding 
velocity of 18 kms ' is used consistent with our non-LTE CO 
modeling of the envelope. For the excitation temperatures, Tex, 
rotational temperatures as computed from Boltzmann diagrams 
are taken (see Table |6]l. Values for the upper energy level and 
line strength (ju^S ) can be found in Table |4] 



4.2.1. Results 

Using the method outlined above, the fractional abundances of 
all molecules (except CO) were determined (see Table [T]) . 

The most uncertain parameters used to derive the fractional 
abundances are Tex, D and x^ (the outer radius). T^x is obtained 
from the rotational diagram analysis, D is taken from the lit- 
erature, and the outer radius of x^, has been adopted differently 
for each individual molecule. We also note that our analysis as- 
sumes optically thin emission, which is not always the case for 
the studied line profiles. The line opacity is expected to be larger 
for higher J rotational transitions, so that lower / rotational tran- 
sitions are expected to better probe the fractional abundance. 
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Table 7. Derived molecular fractional abundance for each transition. 





lid-llftlLlUll 


A hiinrl ixwce^ 




A rM 1 n H Q n 1^ 


V J U. LCI I a\Xi U. a 


I\.Ci Ci CI C 








rpJTiP A^ 


( pa*;p 


( pa*;p 




SiS 


(11-16) 

(20-19) 


1.5x10-" 

U.UA iU 

1.2x10-" 
1.3x10-" 
1.7x10-" 


1 xin'"rm 


3.8x10-' 
1 7vin-' 
3.6x10-' 
4.2x10-' 
5.4x10-' 


Sx10'"rm 




Sin 


("9-1 "1 


1.7x10-5 




6.7x10-" 








(3-2) 


1.5x10-5 




6.2x10-" 




(V\ 




("5-41 


3.8x10"" 




1.5x10"" 








(1-6) 


8.3x10-" 




3.3x10"" 








(8-7") 


1.9x10-5 




7.6x10"" 






sn 




2.0x10-' 




5.2x10-** 










1.1x10-" 




4.3x10-' 








(l-,-(\A 


2.3x10-' 




6.6x10-** 








^"8 '1) 


1.6x10-" 




5.1x10-' 








3^02.) 


1 . / A iU 




J A iU 


S V 1 O^^'rm 
J A iU Clll 






dOiQ-lOnin") 


1.1x10-5 




2.9x10-" 




(T) 




(10oio-9i9) 


1.6x10-5 




5.1x10" 




(2) 






6.0x10-" 




1.7x10-" 








(111 17-16016) 


4.7x10-5 




1.4x10-5 








(43 1-32 2) 


4.8x10-" 




1.5x10-" 








V ^-^2 12 J^-^1 1 1 / 


2.1x10-5 




6.5x10-" 








(■5oT-4n^ 
V-^J 3 ^11} 


2.5x10-" 




8.0x10-' 








n4i m-14j 1 1 ~t 

V lU ^~3 1 1 / 


3.8x10-" 




1.2x10-" 






'"SiO 


(7-6) 
(8-7) 


2.2x10-" 
6.7x10-" 


2xl0'5cm 


8.7x10-' 
2.8x10-" 


5xlO'"'cm 




-'^SiO 


(7-6) 
(8-7) 


6.2x10-" 
1.1x10-5 


2xl0'5cm 


2.5x10-" 
4.6x10-" 


5xlO''cm 




cs 


(2-1) 


4.7x10-' 


lxlO'"cm 


1.1x10-' 


5xlO'"cm 


(3) 




(3-2) 


1.9x10-' 




5.9x10-** 




(1) 




(6-5) 


3.2x10"' 




9.2x10-** 








(7-6) 


2.0x10-' 




6.3x10-** 






HCN 


(1-0) 
(4-3) 


4.9x10-' 
2.3x10-" 


lxlO'"cm 


1.3x10-' 
7.2x10-' 


5xlO'"cm 


(1) 


CN 


N=3-2, J=5/2-3/2 
N=3-2, J=7/2-5/2 


9.8x10-" 
2.3x10-' 


lxlO'"cm 


3.1x10-** 
7.1x10-** 


5xlO'"cm 





Notes. For some molecules, other line transition were searched in literature. 

References. (l) lBuiarrabal et aljjl994h : (2) IOmont et alj ( ll993h : (3) lLindavist et all ( 1198^ . 

Table 8. Comparison of the derived molecular fractional abundances with other published results. 





CS 


HCN 


SiO 


SiS 


SO 


SO2 


CN 


This worlc (case A) 


3.0x10-' 


1.4x10-" 


1.3x10-5 


1.3x10-" 


7.8x10-' 


1.4x10-5 


1.6x10' 


This work (case B) 


8.1x10-** 


4.3x10-' 


5.1x10-" 


3.7x10-' 


2.7x10-' 


4.2x10-" 


5.1x10-** 


(1) 


1.0x10-' 


9.8x10-' 


1.7x10-5 


4.4x10-' 


2.6x10-" 






(2) 


3.0x10-' 


6.0x10-' 




7.0x10-' 








(3) 






3.0x10-" 




1.8x10-" 


4.1x10-" 




(4) 


2.9x10-' 


1.4x10-' 


3.2x10-5 


3.5x10-" 


9.1x10-' 


2.2x10-' 


3 X 10-' 


(5) 


2.8x10-' 


2.1x10-" 


3.8x10-5 


3.8x10-'" 


7.8x10-** 




2.4 X 10-'" 



Notes. In the first part, we list results derived from an observational analysis, in the second part theoretical predictions from chemical models are 
given. Abundances from the chemical models bv iPua ri et al. ( 1999) were selected at a radius of 2.2^, (5 x 10" cm). 

References. (D lBuiarrabal etalJ l ll994h : (2) lLindavist eTZI jl988l) : (3) IOmont et aD ( ll993h : (4) IWillacv & Millar! l l 19971) : (5) lDuari et al1 ( ll999h . 



5. Discussion 

Table [8] and Fig. [14] compare the average abundance of 
each molecule to values found in liter ature. C ompa red to 
observational results from literature (iBuiarrabal et all Il994t 
iLindqvist et al.ll988l:IOmont et al.ll993h . our deduced fractional 
abundances agree within a factor of 3.5 for the smaller outer ra- 



dius (case A) and for the larger outer radius (case B) within a fac- 
tor of 10. Compared to the predicted abundance s fr om theoreti- 
cal chemical models 

bylMlac^AMiiy (Il997h and lDuari et aP 
( Il999h . we found that the predictions are comparable to our de- 
duced values (using the smaller outer radius, case A) within a 
factor of ~3 for CS and SiO. Our deduced value for the SO, 
SiO, CN, and SiS fractional abundances agree with the results of 
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1E-3 




■ This work-A 
# This work-B 

▲ Willacy (1997) 
y Duari(1999) 
4| Willacy (priv. comm.) 



HCN SiO SiS 

Molecular species 

(a) 



HCN SiO SiS 

Molecular species 
(b) 



Figure 14 . Comparison of the m olecular fractional abundances derived in this work (case A and c ase B-study) and values found in 
literature dBuiarrabal et al.ll994tlLindqvist et alJ1988l:IOmont et al.ll993tlWillacv & Millaiill997l;lDuari et aljl999l) . The errors are 
estimated from the abundance variations for the different transitions (see Table[8]l. Panel (a) gives a comparison to other observational 
results, panel (b) to chemical model predictions. 



Table 6. Rotational temperature and beam averaged column den- 
sity. 



Species 


T.O, (K) 


N (cm-^) 


SiS 


85.8(11.1) 


4.46x10'^ (1x10'') 


SiO 


17.1 (1.0) 


8.24x10'" (1x10") 


SO 


27.2 (2.7) 


6.35x10'" (2x10") 


SO2 


67.5 (6.8) 


2.02x10"" (4x10") 


J"SiO 


68.6 (82.3) 


2.48x10'-* (4x10"') 


^^^SiO 


30.0(15.5) 


7.12x10'-* (9x10"*) 


cs 


33.9 (4.7) 


8.89x10'* (2x10"*) 


HCN 


8.3 (0.5) 


2.27x10" (5x10"*) 



Notes. The temperature and the density were determined from the rota- 
tional diagram analysis. The uncertainties are given within parenthesis 
(no systematic errors included). 

["Willacy & Millaii (Il997h . but the predicted values bv lDuari et"ai1 
(Il999h are much lower The SO2 abundance from this work is 
al most two orders of magn itude higher than the value predicted 
bv lWiUacy & Millai] ([T997h . 

As noted above, t he SiS abundance in the chemical models 
of iDuari et al. I ([T999I) is much lower than the observed value. 
The chemical models by iDuari et al ] (Il999h focus on the in- 
ner en velope (within few stellar radii), while IWillacv & Millaii 
0993) studied the chemical processes partaking in the outer en- 
velope. The agreement between our deduced value for the frac- 
tional abundance of SiS and the predictions by I Willacy & Millaii 
(Il997h suggests that SiS is formed in the outer envelope. 

The deduced SO abundance is a factor ~10 higher than 
the inner wind predictions by Duari et al. (1999), b ut agre e 
with the outer wind pr edictions by IWillacv & Millar! (1 19971) . 
IWillacv & MiiiaJ (119971) assumed no SO injection, but only in- 
situ formation. CN is clearly produced in the outer envelope, as 
photo-dissociation product of HCN. 

The abundance of SO2 found bv IWillacv & Millaj (Il997h is 
much lower than the observed ones. A value of 1.4 x 10"^ (case 
A) means that SO2 contains 80% of the solar sulphur value. 



IWillacv & Millaii (1 19971) suggest that SO2 may be formed in a 
different part of the envelope compared to the other sulphur bear- 
ing molecules, for example in shocks in bipolar outflow or in the 
inner envelope. An indication for the typical behavior of SO2 
comes also from the line profiles, e.g. the SO2 (14-14) line is 
clearly narrower and shifted to the red. 

The SiO abundance derived in this study is close to the abun- 
dance predicted by the theoretical chemical models. ICherchneS 
(l2006l) investigated the non-equilibrium chemistry of the inner 
winds of AGB stars and derived an almost constant, high SiO 
abundance ( about 4 x 10~^ before the condens ation of dust). 
IDuari et al.1 ( 1 19991) and IWiUacv & Millad (1 19971) derived 3.8 x 
10"^ and 3.2 X 10"^ for the inner and outer wind, respectively. 
Furthermore, iGonz alez Del gado et al.l (l2003h performed an ex- 
tensive radiative transfer analysis of circumstellar SiO emission 
from a large sample of M-type AGB stars, where they adopted 
the assumption that the gas-phase SiO abundance stays high 
close to the star, and further out the SiO molecular abundance 
fraction decreases due to absorption onto dust grains. Their re- 
sults show that the derived abundances are always below the 
abundances expected from stellar atmosphere equilibrium chem- 
istry. For a mass-loss rate of 4 x lO^^Mo/yTj^ the equilibrium 
chemistry abundance of SiO is -3.5x10"^ Jcherchneff 200^. 
Taking the scenario of depletion due to dust formation into ac- 
count, the higher excitation SiO(8-7) would probe a higher SiO 
abundance. As seen in Table |7] the SiO(8-7) indeed probes a 
higher fractional abundance, although not significantly higher 
than the other lines. 

6. Conclusions 

In this work, we present for the (sub)millimeter survey for an 
oxygen-rich evolved AGB star, being IK Tau, in order to study 
the chemical composition in the envelope around the central tar- 
get. 

An extensive non-LTE radiative transfer analysis of circum- 
stellar CO was performed using a model with a power law struc- 
ture in temperature and density and a constant expansion. The 
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observed line profiles of ^^CO{3 - 2), '^€0(3 - 2), ^^CO{4 - 3), 
and '~CO(7 - 6) are fit very well by our model, yielding a mass- 
loss rate of 4.7 x 10"^ Mo/yr. The line shapes and intensities for 
all '"CO transitions are not much influenced by variations of the 
inner radius, which is understandable since the bulk of the '^CO 
emission is produced in the outer envelope. The intensities for 
the higher excitation CO lines depend strongly on the assumed 
temperature but not on the value of the outer radius. 

For 7 other molecules (SiO, SiS, HCN, CS, CN, SO, and 
SO2) a fractional abundance study based on the assumption of 
LTE is performed. A full non-LTE analysis of all molecules is 
out of the scope of t his observational p aper, but will be pre- 
sented in a next paper ( iDecin et alj|2010l) . This study shows that 
IK Tau is a good laboratory to study the conditions in circum- 
stellar envelopes around oxygen-rich stars with submillimeter- 
wavelength molecular lines. The improved abundance estimates 
of this study will allow refinements of the chemical models in 
the future. 

Molecular line modeling predicts the abundance of each 
molecule as a function of radial distance from the star, although 
some ambiguity about an inner or outer wind formation process 
often exists. To get a clear picture on the different chemistry pro- 
cesses partaking in the different parts in the envelope, mapping 
observations for molecules other than CO should be performed. 
Since most of the submillimeter emission from molecules less 
abundant than CO probably arises from the inner part of the en- 
velope at 2 - 4" meaningful observations require interferometers 
such as the future Atacama Large Millimeter Array (ALMA). 
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